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Abstract. We consider a one dimensional Levy bridge xg of length n and index < a < 2, i.e. a 
Levy random walk constrained to start and end at the origin after n time steps, ig(0) = ig(n) = 0. 
We compute the distribution Pg(A,n) of the area A = X]m=i :c s( m ) under such a Levy bridge 
and show that, for large n, it has the scaling form Pg(A, ri) ~ n~ 1 ~ 1 / a F a (A/n 1+1 / a ), with the 
asymptotic behavior F a (Y) ~ Y~ 2 ( 1 + a ) for large Y. For a = 1, we obtain an explicit expression 
of Fi(y) in terms of elementary functions. We also compute the average profile (ia(m)) at time 
m of a Levy bridge with fixed area A. For large n and large m and A, one finds the scaling form 
(isH) = (m/n,A/n 1+1 /°'), where at variance with Brownian bridge, H a (X, Y) is a non 

trivial function of the rescaled time m/n and rescaled area Y = A/n 1+1 / a . Our analytical results 
are verified by numerical simulations. 
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1. Introduction 

Random walks, and the associated continuous-time Brownian motion (BM), are ubiquitous in nature. 
As such, they are not only the cornerstones of statistical physics [TJ HJ [3] but have also found many 
applications in a variety of areas such as biology @] , computer science E] and finance [7] . Continuous 
time Brownian motion is simply defined by the equation of motion 

x(0) = 0, ^ = »?(<), (1) 

where r)(t) is a Gaussian white noise of zero mean (r](t)} = and short range correlations {ri{t)r]{t')) = 
2D5(t — t') where D is the diffusion constant (in the following we set D = 1). An interesting variant 
of Brownian motion in a given time interval [0,T] is the so called Brownian bridge xsit) which is a 
Brownian motion conditioned to start and end at zero, i.e. xb{T) — x b {0) =0. Here we focus on two 
interesting observables associated with this bridge, namely 

• the distribution of the area A under the bridge (see Fig. [l]a)) 

A = [ x B (t)dt, (2) 



which is obviously a random variable, being the sum of (strongly correlated) random variables. 
For the Brownian bridge, the distribution of A can easily be computed using the fact that x B (t) 
is a Gaussian random variable. This can be seen from the well known identity in law [2J 

x B (t):=x(t)-^x(T), (3) 

where x(t) is a standard Brownian motion ([!]). For the Brownian bridge, A is thus also a 
centered Gaussian random variable. A direct computation of the second moment (A 2 ) yields 
straightforwardly 



^,T)^-^exp(-^). (4) 

the average shape of the bridge, (ifi(t)) for a fixed area A (see Fig. [l]b)). For the Brownian 
bridge, it takes a simple form [5] 

(x B (t)) = if(l) , f(x) = 6x(l-x). (5) 



T V T , 

The distribution of the area under a Brownian bridge Q is a standard result and its extension to 
various constrained Brownian motions has recently attracted much attention EH EQ E3 EH EH E3 • 
For instance, the distribution of the area under a Brownian excursion (i.e. a Brownian motion 
conditioned to start and end at and constrained to stay positive in-between), the so called Airy- 
distribution, describes the statistics of the maximal relative height of one-dimensional elastic interfaces 
[9] EQ [15] . Another example is the area A under a Brownian motion till its first-passage time t f [10] , 
which has an interesting application to the description of the avalanches in the directed Abelian 
sandpile model proposed in Ref. [TBJ, such that tt relates to the avalanche duration and A to the size 
of the avalanche cluster. Related quantities were recently studied in the statistics of avalanches near 
the depinning transition of elastic manifolds in random media |17j . On the other hand, the average 
shape of random walk bridges with a fixed area A ([5| has been studied some time ago in the context 
of wetting [8] to prove the validity of the Wulff construction in 1 + 1 dimensions and more recently in 
the context of mass transport models [18] . In these models, where the transport rules depend on the 
environment of the departure site, the steady state has a pair-factorized form 19J, which generalizes 
the factorized steady states found in simpler system like the zero range process [5D] [STJ [22] . As the mass 
density crosses some critical value, the system exhibits a condensation transition which is governed 
by interactions, which in turn give rise to a spatially extended condensate. It was shown in Ref. 18J 
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Figure 1. a): The blue line is the trajectory of a Brownian bridge xg(t) on [0, T], xg(0) = xg(T) = 
and in red is the area A under this Brownian bridge, b): The typical trajectory of a Brownian 
bridge xg(t) (in blue) with a fixed area A (in red) on the interval [0, T\. The solid black line is the 
average profile (ifl(t)) as given by Eq. j5|. 

that the shape of this condensate can be described by the average shape of a random walk bridge with 
fixed area, and the results of Ref. [5] were recovered. 

While these quantities are well understood for a Brownian bridge, much less is known for the case 
of a Levy bridge. The aim of the present paper is to compute the distribution of the area and the 
average shape for a fixed area A in that case. To this purpose, it is convenient to consider a random 
walk x(m), in discrete time (see Fig. [2] a)), starting at x(0) — xq at time and evolving according to 

x(m) — x(m — 1) + 77(771) , (6) 

where 77(771) are independent and identically distributed (i.i.d.) random variables distributed according 
to a common distribution cf>(r)). Here we focus on the case where 0(77) = S a (r]) where S a (n) is a 
symmetric a-stable (Levy) distribution. Its characteristic function is given by J^' 00 S a (r])e' lkri dr] = 

e _ ' fc ' . In particular, for large 77, the distribution of 77 has a power law tail 4>(r}) ~ n~( 1+a \ with 
< a < 2. The Levy bridge ig(m), on the interval [0, n], is a Levy random walk conditioned to start 
and end at 0, i.e. #3(71) = xb(0) = 0. In the following we will compute the distribution Pe(A,n) of 
the area A under the bridge of length n, i.e. A = J2m=o XB ( m ) an< ^ ^ ne avera g e shape of a bridge 
(ig(ro)), with < m < n, for fixed A. Here we consider the natural scaling limit where xb ~ n 1 /", 
while A ~ 77 1+1 / a , whereas the aforementioned previous works [8j [T8] focused on a different scaling 
limit which, for Brownian motion, corresponds to xg ~ y/n and A ~ n 2 . Note that, even in this 
natural scaling limit, the identity in law valid for the Brownian bridge ^ does not hold for a Levy 
bridge [121 HI] and one thus expects the distribution of the area A to be non trivial. Our results can 
be summarized as follows: 



• for a free Levy random walk starting at Xq — 0, one finds that the distribution of the area P(A, n) 
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P(A,n) = —S a (— 

In \7« 



In 



1/a 



(7) 



where S a (x) is a symmetric a-stable distribution (21 ). For a Levy bridge, one finds, in the scaling 
limit A —> oo, n — > oo, keeping A/n 1+1 / Q fixed, that the distribution of the area Pb(A, n) takes 
the scaling form 



(8) 



where F a (y) is a monotonically decreasing function, with asymptotic behaviors 




F a (Q) , Y^O, 



(9) 



where -F Q (0), see Eq. (36 1, and a a , see Eq. (37 1, are computable constants. For a = 1, we obtain 
an explicit expression for -Fi(y) in terms of elementary functions (C.6). 

on the other hand, in the aforementioned scaling limit, one obtains the average profile (is(m)) 
for a Levy bridge as well as the average profile (x(m)) for a free Levy walk with a fixed area A. 
For (x(m)} one obtains a farely simple expression 



(x(m)) 



a + 1 








1 - 




a 







(10) 



For a Levy bridge, the expression is more involved. For generic a one finds the scaling form 



A 



(11) 



which, in general, has a non-trivial dependence on A. One recovers a linear dependence in A (as 
for the Brownian bridge in Eq. ([5])) only in the limits A — > (68) and ^4 — > oo (69). 



The paper is organized as follows. In section 2, we compute the joint distribution of the position 
and the area under a Levy random walk. In section 3, we use these results to compute the distribution 
of the area A under a Levy bridge of size n while in section 4, we use them to compute the average 
profile of a Levy walk with a fixed area. Finally, in section 5 we present a numerical method, based 
on a Monte-Carlo algorithm, to compute numerically Pb(A,u) and (^s(m)) before we conclude in 
section 6. Some technical (and useful) details have been left in Appendices A,B and C. 



2. Free Levy walk : joint distribution of the position and the area 

We start with the computation of the joint distribution P(x, A, m\xo, Xq, 0) of the position and the 
area after m steps given that x(0) = Xq (see Fig. [2] a)). If we denote by A{m) the area under the 
random walk after m time steps, this random variable evolves according to the equation 

A(0)=x , (12) 
A(m) = A(m - 1) + x(m) . (13) 
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Therefore P(x, A, m\xo, xq, 0) satisfies the following recursion relation: 

P(X, A, 0\XQ, Xq, 0) = S(x — Xq)6(A — Xq) , 
/>oo 

P(x, A, m\xo, Xq, 0) = / P(x — 77, A — x, to — 1\xq, Xq, 0)(/)(ri)dri . (14) 



Introducing (f>(k) = <p(r])e lkr, dr] the Fourier transform of 4>(rf), and thus <p{k) = e ' fe ' for a 

Levy random walk, and Piki, k 2 , m\xo, xq, 0) the double Fourier transform of P(x, A, m\xo, Xq, 0) with 
respect to both x and A, i.e. P(k\, k 2 , m\x , x a , 0) = dx dAP(x, A, m\x , x , Q) e lk i x + lk 2A ^h e 
recursion relation ( |14[ ) reads 

P(k 1 ,k 2 ,0\x ,x Q ,0)^e i ^ +k ^ Xo , (15) 

P(k 1 ,k 2 , m\xo, xq, 0) = 4>{ki + k 2 )P(k 1 + k 2 , k 2 , m — l\x , x , 0) , (16) 

which can be solved, yielding 

n 

P(k 1 ,k 2 ,n\x ,x ,0) = J] $(h + mk 2 )e^ +( - n+1 ^ x ° . (17) 

m=l 

Hence for a Levy walk of index a one has simply 

P(x,A,n\x ,x ,0) = I" ^ [°° * e -KUilW*r e ^M.) e -.MA-(„ + i) IO ) _ (lg) 

J-00 27T J_ 00 2n 



Note that this expression (18) can also be obtained directly by noticing that x(n) — Xq + J>2m=i ^(*) 
and thus A(n) - (n + l)x = IX =i x(m) = Y,m=i YI1L1 ^(0 = ELi mn(n + 1 - to) such that 



poo I \ 

P(x,A,n\x o ,x ,0)=Y\_ / drj(m) <j> [r)(m)} S I x - x - ^ rj(m) J 
x <5 I A - (n + l)x - ^2 m7 l( n + 1 - m ) I ■ (19) 



m— 1 



After a double Fourier transform with respect to x and A, this Eq. ( 19 ) yields immediately the 
expression P(fci, k 2 , n\xo, xq, 0) in Eq. (17). Of course the marginal distribution of the position 
P(x,n\xo,0) and of the area P(A,n\xo,0) are also stable laws. Indeed one has 

P(x,n\x ,0) = j P(x,A,n\x o ,x o ,0)dA= -^S a (j^f) , (20) 
P(A,n\x o ,0)= [ P(x,A,n\x o ,x o ,0)dx = — S a ' S ~ "' ~ ' 



7ri 

where 



in \ Tn 

' " \ 1/Q „1+1/Q 

/ m " ~ 7 7TT7 - , n > 1 , 

^ / (a + 1) 1 /"' 



5„(, )=/ . (21) 



For example, 5i(x) is the Cauchy distribution while S 2 {x) is a Gaussian distribution : 

^) = ~^ 2{ x) = ^=e-4. (22 ) 
Note also the explicit expression 

s a (o)=f e-'^ro+i^ (23) 

./-co 27r 77 

which will be useful in the following. 
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Fig ure 2. a): The area A under a random walk in discrete time, b): Illustration of the method, see 
Eq. |46| , to compute the probability P(x, m\A, n) that the position of the random walker, starting in 
x(0) = is x after m time steps, given that the area, after n time steps, is fixed to A. The light area 
corresponds to A\ , the area under the walk over the interval [0, m] while the shaded area corresponds 
to the area A — A\, the area under the walk over the interval [m,n]. 



We now want to study P(x, A, n\xo, xq, 0) in the limit of large n. The marginal dist ributions in 
Eq. (20) suggest the scaling x ~ n 1 /" and A ~ n 1 ' a+1 . From the expression in Eq. (18) one checks 
explicitly that in the limit n — > oo, keeping X = x/n x t a and Y — Ajn x l aJrl fixed, the joint distribution 
takes the scaling form 



P(x, A,n\x o ,x o ,0) 



1 



n 2/a+l 

where the function G(X, Y\Xq) is given by 



G 



A 



A/a ' n l/ a +l 



Xq 



,1/a 



G(X,Y\X ) = 



dfci 



dk- 



2tt 



?. p -Io \ki+k 2 z\°'d z -ik 1 (X-X„)-ik 2 (Y-X ) 



After the change of variable k2 = k and k\ — kr, we obtain 
r °° dr f 00 dk, 



G(X,Y\X ) 



2tt 



2tt 



\k\e 



-\k\ a -y(r)-ikr(X-Xo)-ik{Y-X ) 



where the function j(r) is explicitly given by 

1 

a -+ 



7(r) 



For a = 2, one has j(r) — 



1 ((-r) a+1 -(-l-rr +1 ) ,r<-l, 

((r + l) a+1 + (-r) a+1 ) , -1< r < 

— — ((r + 1) Q+1 - r a+1 ) , r > . 
a + 1 v ' 

r 2 + r + 1/3 and one finds 
V3 



\r + z\ a dz 



G(X,Y\X ) = ^ exp [-3(Y-X Q )(Y-X)-(X-X ) 2 ] , 
which yields back the propagator of the so-called random acceleration process 



(24) 



(25) 



(26) 



(27) 



(28) 



3. Levy bridge 

In the absence of any constraint for the walker, the area under a Levy random walk is the sum of Levy 
random variables and is thus again a Levy random variable. However, if one considers constrained 
Levy walks, this is not true anymore and the area may become different from a simple Levy random 
variable. In the first subsection, we compute the distribution of the position for a Levy bridge while 
the second subsection is devoted to the distribution of the area under this bridge. 
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3.1. Distribution of the position 

Here we study the Levy bridge {xB(m)}o< m <c n which starts at at time 0, xb(0) = 0, and is 
constrained to come back to after n time steps, i.e. xsin) — 0. In that case, one can compute 
the distribution of the position Ps(a;,m) after m time steps for such a bridge as P B (x, m) = 
P(x, m|0, 0)P(x, n — m|0, 0) /P(0, n\Q, 0), such that in the scaling limit one has 



7T 1 

r(l + a-!) (r(l - r)y/°<~ a VrV« ) ~ a \(1 - t) 1 / c 



Gb(X,t) = —^——±— Sa (X) S a L * „ ) - 



where we have used iS Q (0) = r(l + a _1 )/7r, see Eq. (23). For a — 2 it is easy to see from Eq. (29) that 
xb(iti) is a Gaussian variable. However, for a < 2, the Levy bridge is not any more a Levy random 
variable. For instance, for a = 1 one obtains a non trivial distribution 

GB{X ^-l {T 2 +X 2 ) ^ 1 _ T l ) 2 +X 2 ) ■ (3°) 

It is also easy to see that, for any a < 2 one has the asymptotic behavior 

G B (X, t) ~ c' a T(l - t)X- 2 ^ , X » 1 , (31) 

where is independent of r, which implies that ((xB{m)) 2 ) is well defined for a > 1/2 (of course 
(xfl(m)) = by symmetry for all a). A straightforward calculation shows that 

(x B (m)) _ m r m\ _ aF(2 — a -1 ) 
nV« = Qa n \ n) ' = T(l + a^) 

It is interesting to notice that (x 2 B (m)) /n 2 / a depends on a only through the amplitude a a but the 
parabolic shape in (m/n)(l — m/n) holds for all values of2>a>l/2. Besides a a is diverging for 
a — > l/2 + while one has 5(1) = 1 and 5(2) = 2. One finds, curiously, that it reaches a minimum for 
a non-trivial value a* = 0.74122 . . . for which a(a*) — 0.85264 .... In view of these properties ( 29 30 1 
one expects that, for a < 2, the area under such a Levy bridge has a non trivial distribution, which 
we now focus on. 



5 Q — 1 , a a = -—— — -j— . (32) 



3.2. Distribution of the area 



In this subsection, we consider a Levy bridge, i.e. a Levy walk which starts at the origin Xq = and 
is conditioned to come back to the origin after n steps and we ask : what is the distribution Ps(A,n) 
of the area A under this Levy bridge ? One can obtain Pb(A, n) from Eqs (18 20 ) as 

P(0,An|0,0,0) 



P B (A,n) 



P(x = 0,n|0,0) 



(33) 



Therefore, using P(x = 0,n|0,0) = n _1 / a s S a (0) = n^Z^l + a" 1 )/^, see Eq. g, together with the 
scaling form (24 26) one obtains, in the limit n — > oo, keeping A/n 1+1 / a fixed: 

A 



1 
1 

2r(l + a- 1 ) 



P B(A 1 n) = -^F a 
F a (Y) 



1+1/ c 



dr 



dk 
2tt 



\k\e 



-\k\ a -f(r)-iYk 



(34) 



Using the explicit expression of j(r) above (27), one computes the Fourier transform F a (k) of F a (Y) 
as 



F a (k) 



F a (Y)e lkY dY 



\k\ 

r(i + £) 



(35) 



+ l_r- c < + 1 



dr 



1/2 



-i^[(l/2+r)" + 1 + (l/2-r)° + i ; 



dr 
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For generic a, it seems quite difficult to perform explicitly the integrals over r and k in the expression 



for the distribution F a (Y) in Eq. (34). One can however extract from this expression the asymptotic 
behaviors both for Y — > and Y — > oo. 

Asymptotic behavior for small argument. For small argument, it is straightforward to see 
on the expression ( |34"| ) above that the leading behavior of F a (Y) when Y — > is given by 

F a (Y) ~ F a (0) , Y -> ,with 



F a (0) 




dr 



1/2 



dr 



(( r + l) Q +i_ r «+i)f J Q ((l/2 + r) a + 1 + (1/2 -r)^ 1 ); 



A study of this function ^(0) shows that it is a decreasing function of a on the interval ]0, 2], which 
is diverging when a — > 0. For a = 1 and a = 2, -F a (0) assumes simple values 



Fi(0) = 1 + - = 2.27324... , F 2 (0) = 

7T 

Asymptotic behavior for large argument. 



- = 0.977205. 

7T 



(36) 



The analysis of the large argument behavior of 



F a (Y) is more involved. A careful analysis, left in Appendix B shows that for large Y one has: 

W~^^,y»i, (37) 

/7rr(2 + 2a) tan(a7r/2) 



y-2(H 

2- 2 ( 2 +"). 



r(2 + a- 1 )r(i-a)r(| + a) 



(38) 



When a — > 2, one has from (37) 



Y- (a — 2) . From Eq. ( 37 ) one obtains that the secon d moment 



of the distribution (Y 2 ) is denned only for a > 1/2 where it takes the value (see Eq. (B.27)) 



(Y 2 ) = 



oT(2 



L ) 



12 12r(l + a- 1 ) 



a > 1/2 



(39) 



where the amplitude a a appears in the expression for {x B (m)) computed above (32). This power law 
tail of the area distribution (37) with an exponent 2(1 + a) is quite interesting. Indeed, the area itself 



is the sum of non- identical and strongly correlated variables x b (m) all having a similar power law tail 
also with exponent 2(1 + a) (31). For a > 1/2, their variance is finite and the non-Gaussianity of A 



can a priori be due both to the correlations between the xs(m)'s and to the fact that A is the sum of 
non-identical random variables. To test which of these features is responsible for the non-Gaussianity 
of A, we study the sum of n random variables X m which are independent and such that X m has the 
same distribution as XB{m). Defining S n = J2m=i -^-m an d S 2 t = X)m=i(^m)> ^ ^ s known that S n /Ti n 
converges to a centered Gaussian variable of unit variance if the following condition (known as the 
Lindeberg's condition) is satisfied 

lim — r / x 2 Proba(A m = 2;)d2; = , Ve > . (40) 

n ^°° £2 J\ x \ >fSn 

Intuitively, this Lindeberg condition (40) ensures that the probability that any term X m will be of 



the same order of magnitude as the sum S n must tend to zero (see Appendix A for an example of 
non-identical independent random variables which do not satisfy the Lindeberg condition). In the 
present case of the Levy bridge, one can check (see Appendix A) that if X m is distributed like XB{m) 
then the above Lindeberg condition (40) is satisfied (see Appendix A). Therefore the deviations from 
Gaussianity (37) are purely due to the strong correlations between the positions of the walker 0:5(771.) 's. 

The special case a = 1. In the Cauchy case, a = 1, the integral over k can be done in Eq. (26) 
to obtain 

7 2 (r) - (rX + Y) 2 dr 



G(X,Y\0)= f 

J — t 



(f(r) + (rl + y) 2 )' 



2vr 2 



(41) 
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where the function j(r) (27) takes here a rather simple form 

1 

-"-2 



r < — 1 



7 (r) = < 



1 



, -1 < r < 



(42) 



1 



,r > 



The distribution of the area under a Cauchy bridge is thus given by 



P B (A,n) = -.Fx 

ip ^ 

Fi(Y)= 1 



1 

n 2 
1 2 



.4 



7T 1 + 4Y 2 



V 2 {u 2 + l/4) 2 -Y 2 
{{u 2 + 1/4) 2 +Y 2 ) 2 



du 



Under this form ( 43 ) , one can easily obtain the asymptotic behaviors as 

4 

1 + - , r^o 



Fi(Y) 



1 



(43) 



(44) 



207rF 4 



Y 



oo 



in agreement with the asymptotic behaviors obtained above (|36 37|). In fact the integral over u in the 

Ion 
2 

Y 2 

4 



expression above ( 43 ) can be done explicitly yielding the expression 

1 



Fi(Y) = 



1 



7T 1 + 4Y 2 

2(1 - 8Y 2 ) 



(45) 



7T V (1 + 4F 2 )(1 + 16Y 2 ) (l + IQY 2 ) 



-Re 



(l-4iy)iarctan^(l + 4iy)"^ 



where Re(z) denotes the real part of the complex number z. In Appendix C we show how this 



expression (45) can be written explicitly in terms of elementary functions (C.4 C.5 C.6). 



4. Average profile for fixed area A. 

4-1. The case of a free Levy walk 

We first consider the case of a free Levy walk of constrained area. We compute the probability 
P(x, m\A, n) that the position of the random walker, starting at x(0) = xq — at time 0, is x after m 
time steps given that the area, after n time steps, is fixed to A. From this probability, one obtains the 
average profile as (i(m)) = xP(x, m\A, n)dx. To compute this probability P(x, m\A, n), we divide 
the interval [0, n] into two intervals [0, m] and [m, n]. Over [0, m] the process starts at xq = with area 
Aq = and reaches to x with area A\ (see the light area on Fig. [2]b)). Over the interval [m,n], the 
process starts in x and reaches to if with area A — A\ (see the shaded area on Fig. [2]b)). Therefore 
this probability P(x, m\A, n) can be simply expressed in terms of the propagator P(x, A, n\xo, xq, 0) 
computed above (18) as (see Fig. [2]b)): 

P(x,m\A,n) = -=- r -. — r/ dx F / dAxPix, A 1 , m|0, 0, 0)P(x F , A - A 1: n - mix, x, 0) , (46) 
P{A,n) J^^ 



where we have used the Markov property of the Levy random walk. In the above expression (46), xf 
is the end point of the walk (see Fig. Nb)), which is free here. Hence P(x,m\A,n) is obtained by 
integration over this end point xf ■ Notice that it is normalized according to f_ ao P(x, m\A, n)dx = 1 



(and therefore we have divided by P(A, n) in the expression above (46 1 because the measure is restricted 



to random walks of fixed area A after n time steps). Using the explicit expressions computed above 



(18) one obtains after integration over xf and A\ 
1 r°° dfci 
P(A,n\0,0) J_ 00 ~to 



P{x, m\A, n) 



dh 
2tt 



5ZT=o \ k i+ ly k 2 \ c " -\k 2 \ c " J2" = (T \ v \ a p— ikix— ik 2 (A— (n— m)x) 



(47) 
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In the large n limit, keeping X — x/n 1 / 01 , Y = A/n 1+1 / a and r = m/n fixed one has 
P(x,m\A,n) = -^G(X,t\Y) , 



G(X,t\Y) 



(l + a)~ 1/a 



S a ((a + iy/«Y) 7_oo 2n J _ x 2ir 



dr 



dk 



\k\e 



-|fe| a (7(r,T)+7(0,l— t)) -ikrX-ik(Y -(I-t)X) 



where we have used the expression of P{A, n|0, 0) given in Eq. ( 20 ) and we have introduced 

7(r, t) = \r + z\ a dz , 
Jo 

which is a generalization of the function j(r) = 7(r, 1) in ( [26] ) . It reads 

r < —t , 

-t < r < , 



(48) 



(49) 



7(r,r) = < 



1 


(X 


+ 


1 




1 




a 


+ 


1 




1 




a 


+ 


1 



I ((r + r)«+ 1 + (-r) Q+1 ) 



(50) 



\a + l _ r a+l\ 



j , r>0. 

We can now compute (x(m)) which, in the large n limit, takes the scaling form 

(i(m)) 



n X/a 

K{t,Y) = 



m A 
~> n i+i/ a 

(1 + a)-^ a 
S a ((a + iy/*Y) 



(51) 



dXX 



dr f°° dk 



2tt 



2tt 



\U e -\k\ a g(r,T) e -ikrX-ik(Y-(l-T)X) 



with g(r, r) = 7(7", r) + 7(0, 1 — r). This function h a (T, Y) can be written as 



K{t,Y) 



(1 +a)- 1 / a 



dX 



dr f°° dk \k\ 



-\k\ a g(r 



-ifcrX-ifc(Y-(l-r)X) 



(52) 



5 Q ((a + 1)1/^)7^ J_ 00 2n J„ 00 2n(-ik) 

which suggests to perform an integration by part in the integral over r, yielding (one can check that 
the boundary terms vanish) 



K{t,Y) 



(l + a)- 1 /" 
S a {{a + iy/*Y) 



dX 



dr r°° dk \k\ 1+a dg{r,r) 
2^ J_ 00 2^ (-ik) dr~ 



e ~\k\ a g(r,r) £ -ikrX-ik(Y -(1-t)X) 



(53) 



On this expression (531, the integral over X can be done yielding simply a delta function of r, namely 
2ir\k\~ 1 8{r — (1 — tJ). This allows us to perform then the integral over r to obtain 



h(rY)- ( 1 + Q )' 1/a j (Mill 



dk 
2^ 



(54) 



r=l— r ^ -00 

where we have used the relation 7(1 — r, r) + 7(0, 1 — r) = 1/(1 + a). It is then easy to check that 



27r a L 



so that finally one obtains the simple result 

Mr,r)=r — (l-(l-rn 
a 



(55) 



(56) 



where we have used ( da ^ T ^ 
function of r. 



1 — (1 — r) Q . In Fig. (3 1 a), we show a plot of h a (r,Y)/Y as a 
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a) T b) 



Figure 3. Average shape of a free Levy walk in a) and of a Levy bridge in b). a): Plot of h a (t, Y)/Y 
given in Eq. | |56| l as a function of r for a = 2 (solid line), a = 1 (dotted line) and a = 1/2 (dashed 
line), b): Plot of Hi(t, Y)/Y as a function of r in the limit A — > (solid line), as given by Eq. 
and in the limit A — > oo (dotted line), as given in Eq. j69| l. For comparison, we have also plotted 
the parabola (dashed line) corresponding to the result for a Brownian bridge, i.e. a = 2, as given in 
Eq. |5| 



4-2. The case of a Levy bridge 

Here we consider a Levy bridge, i.e. a Levy random walker starting in at initial time and constrained 
to come back to the origin after n time steps. We compute the probability Pg(x,m\A,n) that the 
position of the random walker is x after m time steps given that the area, after n time steps, is fixed 
to A. From this probability one obtains the average profile as (xB(m)) — xPb(x, m\A, n)dx. 
This probability PB(x,m\A,n) can be expressed, as in Eq. (46) in terms of the propagator 
P(x, A, n\xo, a 

Pb{x, m\A, n) 



P(x, A, n\xo, xq, 0) computed above in Eq. (18) as 
1 

P(0,A,n\0,0,0) 



P(x, A x ,m\Q, 0, 0)P(x, A-A x ,n- m|0, 0, 0)cL4i 



(57) 



where we have used the Markov property of the Levy random walk. It is normalized according to 
J_ PB(x,m\A, n)dx = 1 (and therefore we have divided by P(0, A, n|0, 0, 0) because the measure is 
restricted to bridges of fixed area A). Using the explicit expressions obtained above (18) one has 

r0 ° dfci 



Pb(x, m\A, n) - 

In the large n limit, keeping X 
Pb(x, m\A, n) - 

G b (X,t\A) = 



1 



P(0,A,n\0,0,0)J-oo 2^ 



27 



dh 
2^r 



2 ^—i{h\-{-k'i)x—ikA 



Yl™=i \ki+vk 2 \ a ~Yl™=r \k[+vk2\ c 



X e 

x/n 1 / 01 , Y — A/n 1+1 / a and r = m/n fixed one has 
1 



(58) 



l/a 



n 



G b (X,t\Y) 



x e 



r(l + a-i)F a (F) 

-|/c| ct [7(r,r)+7(r',l — t)] 



00 dr 
-oo 2^ 



00 dr^ 
-oo 2^ 



dk 
2tt' 



^2 e -ik{r+r')X-ikA 



We can now compute (is(m)), which in the large n limit takes the scaling form 



where 



(x B {m)) 

n l/a 



H a {r,Y) 



= Ha 



A 



,1+l/c 



dk 



2Y{l + a )F a (Y) J_ x 2tt 



dr|fc| Q - 1 /fcc) r 7(r,T)e- |fc|Qf ^ r) - l ' £y 



(59) 
(60) 

(61) 
(62) 



Area distribution and the average shape of a Levy bridge 

1 d 



12 



2T(l + a)F a {Y) dY 



dk 
27 



dr|fcr- 1 a i .7(r,T)e- |fc|af(r ' r) - lfeY 



where we have introduced the notation T(r, r) = 7(7*, r) + 7(— r, 1 — t) which we compute 
straightforwardly from Eq. (50) as: 



r(r, r) = 7(r, r) + 7(— r, 1 — r) = < 



1 


a 


+ 1 




1 


a 


+ 1 




1 


a 


+ 1 



r + 1 - t) q+1 - (-r - r) a+1 ) , r < -r 

- ((r + t) q+1 + (-r + 1 - t) q+1 ) , -r < r < 1 - r 

— — ((r + t)" +1 - (r - 1 + t) q+1 ) , r > 1 - t . 
^ a + 1 v y 



(63) 



Note that this function T(r, t) satisfies the identity 

f(l -T-r,r)= f \z - r\ a dz = j(-r) 
Jo 

independently of r. 



(64) 



For generic a, the expression above (62) is quite difficult to handle. For a — 2 (Brownian motion) 
and a = 1, further analytical progress is however possible. For a = 2, one has <9 r 7(r, t) = r(2r + r) 
and r(r, r) = r 2 + r(2r — 1) — r(l — t) + 1/3 and therefore one checks 



d r j(r, r)e- fc2r(r ' r) dr = ^Je"^T(l - t) 



(65) 



Using this identity (65) and integrating over r in the expression above (62), and using F 2 {Y) = 
-\/3/7rexp (— 3Y 3 ) one obtains for a = 2: 

# 2 (r,Y) = 6Yt(1-t) . (66) 

Another inter estin g case where analyti cal progress is possible is a = 1. Given the expres sion 
of 7(r, t) in Eq. (50) and T(r, t) in Eq. (63), one observes that the integral over r in Eq. (62) 
gives rise to 4 different terms, corresponding to r 6] — 00, — r], r € [— r, 0], r G [0, 1 — r] and finally 
r 6 [1 — 7", +00 [. For a = 1 it turns out that the first and fourth terms, corresponding to r g] — 00, — r] 
and r G [1 — r, +00 [ do cancel each other (which is the case only for a = 1) resulting in the following 
expression: 



#i(t,Y) = 



Y 



nT(l + a)Fi{Y) 



(2r + rW 



T(r,r) 



(r(r,r) 2 + F 2 ) 2 



-dr + t 



r(r,r) 



(r(r,r)2 + y2)2 



-dr 



with f (r, t) = i((r + r) 2 + (r + r — l) 2 ). In the asymptotic limit Y — > one obtains 



Hx{t,Y) 



Y 



4 + 7T 

In the opposite limit Y — > 00 one obtains 

.10 



3tt - 2 H f + 12(2t - 1) arctan (1 - 2r) 

1 + 2t(t — 1) 



Hi(t, F) ~ Y-t(1 - t)(2 - r(l - r)) 



(67) 



(68) 



(69) 



Note that although the two functions of r entering these asymptotic expansions in Eq. ( |68| and 
Eq. (69) have very different analytical expressions, they are actually quite close to each other on the 
interval [0, 1] (see Fig. ^ b)). 



5. Numerical results 



We now come to numerical simulations of Levy bridges. As mentioned above, one can not use 
the relation above ([3]), which is only valid for a = 2 [23] to simulate a Levy bridge. Instead, we 
consider the joint probability distribution function (pdf) of the increments rj(m) for a Levy bridge 
of size n. Indeed, these increments are independent random variables, distributed according to 4>{rj) 
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with the global constraint that x(n) = J2m=i vi 71 ) — 0- Therefore the joint pdf of the increments 
p B (r?(l),7/(2), • • ■,ri(n)) is simply given by 

n / n \ 

P B (77(13,77(2), • • • , 77(77)) x j { o [77(777)] <5 ]T 77(777) (70) 



m— 1 



oc exp 



.m— 1 




This joint distribution can thus be considered as a Boltzmann weight with an effective energy 
E = — X)m=i m [0 W 771 )]] an d effective inverse temperature (3 — 1- This thus leads us to use a Monte- 
Carlo algorithm, with a global constraint, to generate "configurations" of the increments distributed 
according to the distribution above ( |70"| . We implement it in the following way. We start with a 
random initial configuration of the 77(771) 's which satisfies the global constraint Y] , 77(777) = (it can 
also be 77(777) = 0, for all to). At each time step we choose randomly two sites i and j among 1, 2, • • • , 77 
and the simple following moves are proposed 

77(1) -> r)'(i) = 77(7) + A77 , 

vti) -> rftf) = 77O') - A77 , (71) 
such that the global constraint of zero sum is automatically satisfied. This move is then accepted, in 
Metropolis algorithm that we use here, with a probability given by 



This Monte Carlo algorithm is thus very similar to the Kawasaki dynamics for ferromagnetic spin 
systems relaxing towards equilibrium with a conserved global magnetization [26 . Once the increments 



?7(fc)'s are generated according to this joint probability (70), we can generate the random walk bridge 
£5(777) = Ysk=i v(k) an d compute the distribution of the area A — J2m=i x B( m ) under the Levy 
bridge. In Fig. [4] a), we show a plot of this distribution P B (A,n) for a — 1 and 77 = 100. To compute 
it we have first run 10 Monte Carlo steps to equilibrate the system and the distribution was then 
computed as an average over 10 7 samples generated in the time interval [10 7 ,2.10 7 ]. In Fig. [ZJ we 



also show a plot of the exact explicit expression for Fi(Y) given in Eq. ( C.6[ ), showing a very good 



agreement with our numerics. We have also computed numerically this distribution for other values 



of a G]0, 2[, showing a good agreement with the power law tail obtained in Eq. (37|. Note however 
that for small a, it is actually quite difficult to equilibrate the system such that a precise estimate of 
the exponent characterizing the power law tail of Ps(A,n) is quite difficult for a < 1. 

We can use a similar Monte Carlo approach to generate a random walk bridge with a fixed area 
A = Y^m=i( n + 1 ~ 777)77(777). In that case, the joint pdf of the increments Pb (77(1), 77(2), • • • , 77(77)) is 
simply given by 

n I n \ n 

P B (77(1), 77(2), 77(77)) cx J]_ [77(777)] 5 ) ^(^(77 + 1-777)77(777) - A) . (74) 

m— 1 \m— 1 / m— 1 

We start with an initial configuration of the r)(m)'s which satisfies the both global constraints. In 
practice, we start with 77(777) = 6A(N + 1 — 2m) /((N — 1)N(N + 1)). Then, to satisfy both constraints 



(74), at each time step we choose randomly three sites i, j and k among 1,2, • • • ,77 and the simple 
following moves are proposed 

77(7) -> 7/(7) = 77(7) + A77 , 

7 — k 

vti) -> v'U) = vU) + 7^r~ Ar ? ' 

T)(k) -¥ 77' (k) = n(k) + 7 — -A77 . (75) 
k-j 
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Y = A/n 



0.2 0.4 0.6 0.( 

r m/n 



Figure 4. a): The squares represent the numerical data for n 2 Pg(A, n) as a function of Y = A/n 2 
for a Levy bri dge o f index a = 1 of length n = 100. The solid line is our exact expression for F\(Y) 
given in Eq. JcTm. b): The squares represent our numerical data for n(xg (m))/A as a function 



of t — m/n for a Levy bridge of index a = 1 of length n = 
corresponds to our asymptotic result for large A given in Eq. 
the result for the Brownian bridge in Eq. |5j. 



100 and A/n 2 = 20. The solid line 
||69| while the dotted line represents 



Note that to converge to the correct probability measure ( 74 ) one has to choose i] either positive or 



negative with equal probability. This move ( 75 ) is then accepted with a probability P^ given by 

<^Wfa'(i)W(*0 



Pijk = min I 1 



min(l, exp (-AE)) , AE = log 



<t>W)]<t>W)]<t>[r){k)] 



(76) 



(77) 



Once the increments tj(A;)'s are generated according to this joint probability (74), we can generate the 
random walk bridge is(rn) = Y^k=i wft) with fixed area A and compute the profile (xB(m)). In Fig. [4] 
b), we show a plot of this average profile for a — 1, n — 100 and A/n 2 ~ 20. To compute it we have 
first run 10 7 Monte Carlo steps to equilibrate the system and the average was then computed over 10 7 
samples generated in the time interval [10 7 ,2.10 7 ]. In Fig. [4] b) , we also plot, with a solid line, our 
asymptotic result in Eq. (69), showing a relatively good agreement with our numerics (note that here 
A/n' 2 — 20). On the same plot, Fig. [4] b) , we also show in dotted line, the result for the Brownian 
bridge which is independent of A. It is quite remarkable that these two profiles are very similar 
which show that the global constraints that we impose here have strong consequences on the statistics 
of the Levy random walk. 



6. Conclusion 

To conclude, we have studied two main properties of a Levy bridge ig(m) of length n : (i) the 
distribution Pb(A, n) of the area under a Levy bridge and (ii) the average profile (xB(rn)} of a Levy 
bridge with fixed area A. 

• For Pb(A, n) we have found the scaling form, valid for large n, Ps(A,n) ~ ?i _1_1 / Q F a (Y) 
with an interesting power law behavior F a (Y) ~ y- 2 ( 1 + Q ). For a = 1, we have obtained an 



explicit expression for F±(Y) in terms of elementary functions (C.6). We have also shown, using 



the Lindeberg condition that the non-Gaussianity of Pb{A, n), for a > 1/2, is due only to the 
correlations between the positions of the walkers xs(m)'s. 
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• For the average profile, (Jb(to)), we have found the scaling form (Jb(to)) ~ 
n l l a H a (m/n, A/n 1+1 l a ) where, at variance with Brownian motion, H a (X, Y) is a non trivial 
function of the rescaled area Y. For a = 1, we have obtained simple analytical expressions for 
Hi(X, Y) in both limits Y — > and Y — > oo. In particular, we have shown that the average profile 
of the Levy random walk with a fixed area is not very far from the profile of a Brownian bridge 
with fixed area. 

• We have finally compared our analytical results with Monte Carlo simulations of these Levy 
random walks with global constraints. 

In view of recent developments in the study of area distributions for variants of Brownian motions 
[T51 [3 [TTJ [TS] , it would be very interesting to extend the results presented here to other constrained 
Levy walk, including in particular Levy random walks conditioned to stay positive (Levy excursion), 
which is a challenging open problem. 
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Appendix A. On the use of the Lindeberg condition 

Let us consider independent and non-identical random variables X\ , • • • , X n which have the same 
distribution as the Levy bridge x b (to) ( 29 1 , i.e. 

7T Tl^/ a / X \ ( X \ 

Proba(X m = x) = r(1 + a -i/2 )m i/ Q(n „ m) i/a 5 « Krr^x) Sa (j„_ m )i/a ) ' ( A ' X ) 

For q > 1/2, ct^j = (X^) is well defined and one has (32 1 

a 2 m = = (4 (m)) = a a n 2 / a - 2 m(n - to) . (A.2) 

Given that the variables Xj are not identical, one can not apply directly the Central Limit Theorem. 
However, one can show that these random variables X m do satisfy the Lindeberg condition which 
guarantees that their sum A n — X)m=o -^m ^ s distributed according a Gaussian distribution in the 
large n limit. Let us first introduce £ 2 

n 

K = £ = y- 2/Q ~ 2 (" + !)"(" - 3) ~ Y nVa+1 ' n >> 1 ' (A - 3) 

m— 1 

which implies £„ ~ tt, 1 / 2 + 1 / Q! for large n. To apply the Lindeberg condition, we need to estimate for 
any e > 

(X 2 n ) e = { x 2 Proba(X m = x)dx (A.4) 

J\x\>eS n 

A/a 



2 ; nl/a r x 2 s ( x \s ( x \ 

L(l + a- 1 /2) TO i/a( n _ m )i/a J a \ m i/*J a \{n~m) 1 / a ) 



dx 



~ c> 2/Q n-( 3 / 2+Q )m(n - to) , (A.5) 
where c' a is independent of to and n. Therefore one has 

Tl=^xl) t ^ n - {a -i/ 2) (A 6) 

For random variables for which the above ratio ( |A.6[ ) goes to zero in the limit n — > oo (which is the 
case here for a < 1/2), a theorem due to Lindeberg (thus called the 'Lindeberg condition') [2J, says 
that their sum A n /Yl n = S^ 1 X)m=o^ m ^ s distributed, in the limit n — > oo, according to a Gaussian 
distribution of unit variance. The fact that, for a Levy bridge, the area is not a Gaussian distribution 
(37) is thus, for a > 1/2, a consequence of the correlations between the random variables Xb{to). 

To conclude this paragraph, we discuss a simple case where the Lindeberg condition does not hold. 
Consider the case where X\ , ■ ■ ■ , X n are independent random variables distributed according to [37] 

Proba(A r „ =x) = me~ mx , (A.7) 
such that one has 

(X m ) = - , ((X m - (X m )) 2 ) = \ . (A.8) 

TO TO 

Then in that case one has immediately 

n n -i 2 

K = Y. d X m - (X m )) 2 ) = £-2^>^oo. {A.9) 

Til O 

m— 1 rn—l 

One computes straightforwardly, for e > 

p-eS n m 

(A m - (A m )) 2 Proba(X m = x)dx = — (l + TO 2 (e£„) 2 ) , (A.10) 

eS„ m 

such that here one has 

-, oo /.oo 



v ., / , (X m - (X m )) 2 Proba(X m = x)dx -> — > — (1 + (e'to) 2 ) > , (A.ll) 

^ m =1^6S n T TO 



with e' = e7r 2 /6. Therefore the Lindeberg condition (40) does not hold here. In fact, it can be shown 
that the distribution of the variable (S„ — X)fc=i k~ 1 )fE n converges to a Gumbel distribution [2"T| . 



Area distribution and the average shape of a Levy bridge 
Appendix B. Asymptotic behavior of F a (Y) for large Y 



17 



To analyse the large argument behavior of F a (Y), we analyse the small k behavior of its Fourier 
transform F a {k) given in the text in Eq. (35 1: 



F a {k) 
F a ,i(k) 

F a , 2 (k) = pTj 



AYF a {Y)e lkY = F a>1 (k) + F Qj2 (fc) , 



\k\ 
k\ 



1/2 



i£[ (H -l)«+*-r«+i] drj 
.^[( 1 /2+r) a + 1 +(l/2-r) D '+ 1 ] dr 



(B.l) 
(B.2) 
(B.3) 



The analysis of the small k behavior of F a ^(k) is simply obtained by expanding the exponential under 
the integral. It yields straightforwardly: 



F a , 2 {k) 



1 



r(i + a - 



\k\' 



\k\ 1+2c 



27rr(2 + a) 



(a + l)(a + 2) 2(a + l) 2 I 3 + 2a r(|+a)2 3 + 2 



■0(\k\ 



(B.4) 

The asymptotic expansion of F ai i(k) is a bit a more subtle. To get the two first terms of the 
expansion, one performs the change of variable z = \k\r and then expand (r+ 1) Q+1 using the binomial 
formula, 

1 <T 

(B.5) 



This yields 



F aA (k) 



1 



Til + a- 1 ) 



H fc l''° 1 dr + 0(fc 1+ ") , 



(B, 



where 77 > is yet unknown (see below). From this expression (B.6 1, one immediately obtains the two 
first terms of the expansion of F ay i(k) as 

2r(TT^j + < t,+, '>' 



F a ,i{k) = 1 



(B.7) 



Combining Eq. (B.4) and Eq. (B.7) one sees that the first non-trivial term , pro portional to |fc| cancel 
in F a (k). Therefore, one needs to develop F a> x(k) beyond the first terms (B.7). To this purpose, we 
need to study separately the cases < a < 1/2, 1/2 < a < 1 and 1 < a < 2. 

Appendix B.l. The case < a < 1/2 

Let us first analyse the term F a i(k) which we decompose as 



F a ,i(k) = B^k) + B 2 (k) + B 3 (k) 



B 1 (k) = 
B 2 (k) = 
B 3 (k) = 



\k\ 


r(i 






\k\ 


r(i 






k 


r(i 


+ a~ 1 ) 



dre -ifcr('- Q +f ifek Q_1 ) _ 



(B.8) 

(B.9) 
(B.10) 



It is easy to expand Bi(k) for small k as 

\k\ 1+a 2-2«+ 2 |fc| 1+2Q 



r(l + a- 1 ) (a + l)(o + 2) 2r(l + a- 1 )J 



dr 



T ((r + 1) 



o+l _ r a+l\ 



+ 0{k l+Aa ) 



(B.ll) 
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To expand B 2 (k), one checks the asymptotic behaviors, for large r 

-^—Hr + l) a+1 - r a+1 ] - (r a + -r^ 1 ) = 0(r a - 2 ) 
a + l LV ' 1 v 2 ' v ' 



1 



- (r a + -r*- 1 ) 2 = 0(r 2a - 2 ) 
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(B.12) 



(B.13) 



so that, for a < 1/2 one can safely expand the exponentials in the integrand of B 2 {k) (B.9) up to 
second order to obtain 



B 2 (k) 



dr 
dr 



1 



T [(r + 1) 
1 



(r a + 2^ _1 ) 



a+l _ r a+l 



r(i + a-i)A 
\k\ 1+2a f°° 

+ 2r(l + a-!) J 1 " \\a + l 

+ (| fc |min(2,l+3a)) _ 

Summing up the contributions from Bi{k) and B 2 {k) and performing the integrals yields 



(r a + 2 rQ_1 ) 2 



(B.14) 

(B.15) 
(B.16) 



B 1 (k)+B 2 (k) = 



\k\ 



l+a 



1 1 



r(l + a- 1 ) \(a+l)(a + 2) a + l 2 



1*1 



1+2q 



1 1 



(B.17) 
r(-3 - 2a)r(l + a 



2r(l + a- 1 ) [2 l + 2a (l + a 2 )(3 + 2a) 8a - 4 

+ c ,(| fc |mm(2,l+3a)) ^ 

which we have carefully checked using Mathematica. 



r(-«) 



(B.18) 



Let us now expand B 3 (k) (B.10) for small k. It is easily seen from Eq. (B.7) that the first terms 

o{\k\ 1+tl ) , 



of this expansion are indeed given by 

1*1 



B 3 (k) = 1 



(B.19) 



2r(l + a" 1 ) 

with /i > 0. To go beyond the lowest orders, we first perform a change of variable x — \k\r and then 
compute B' 3 '(k) and finally expand it for small k. This yields 



B 'i [k) = r(l + a-') { ai \ + (1 + % ))*°"' e " <1+S "° + . 



2 Z 100 



V^" 1 ) da; 



(B.20) 
(B.21) 



where 2(a — 1) < —1 for a < 1/2. One then obtains the small k behavior of ls(*) by simply expanding 
the term e~^ kx in the integrand. This yields, to lowest order 

1 /a\ 2 



h(k] '= : T3to(f) fc2 ^ 1 + °( fc3Q " 1 ) 

on 
1*1 



From Eq. (B.20) and Eq. (B.22I, one obtains straightforwardly 

I Lll+a 



B 3 (k) = 1 - 



1 1 



r(l + a- 1 ) r(l + a- 1 )\ v 2 a + l 
|fc| 1+2Q / l + a/2 (I 



r(l + a-!) V 2(1 + 2a) \2 



+ (l + a/2) + 



1 /a\ 2 



2a(l + 2a) 1 - 2a \2 



(B.22) 

(B.23) 
(B.24) 



Fina lly, combining Eq. (B.17) and Eq. (B.23) together with the small k expansion of F a<2 (k) above 
(B.4), one sees that the term proportional to |/c| 1+Q actually cancels, yielding 

F a {k) = l + c a \k\ 1+2a + 0(k 1 + 3a ) , 

1 2- 2 ( 2 + Q )7r 3 / 2 tan(a7r/2) 



T(l + a- 1 ) cos (a?r)(l + a)r(-a)r(5/2 + a) 



(B.25) 
(B.26) 
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This singular behavior of F a (k) for small k (B.25) yields the power law behavior of F a (Y) for large Y 
W"^. (B.27) 



y2(l+ 

1 . „ . 2- 2 ( 2 + a )0i : r(2 + 2a)tan(a7r/2) 

-r(2 + 2a) cos (a7r)c Q = v — K — '—^ 

ir K ' V ' r (2 + a- 1 r 1 -a )T § + a) 



(B.28) 



Appendix B.2. The case 1/2 < a < 1 

This case can be studied along the same line as above except that in that case, l + 2a > 2 and therefore 
one has to handle with care the analysis of terms which are proportional to k 2 , while the coefficient 
proportional to k 1+2a has the same form (B.25). We will not repeat the analysis and simply give the 
result. One finds that F a (k) behaves for small k as 

b, 



F a {k) 



1 



^k 2 +c a k 1+2a 



0{k 



min (3,l+3a)'i 



aT{2~a~ 1 ) 
12r(l + ar T ) 



(B.29) 



The expression for b a given above (B.29) yields the expression for (Y 2 ) given in the text in Eq. (39). 



Appendix B.3. The case 1 < a < 2 

In this case one can again perform a similar analysis but in this case one has 1 + 2a > 3. And therefore 
one has to handle carefully the term proportional to \k\ 3 . A quite lengthy calculation shows that this 
term actually vanishes for a > 1, while the coefficients of the terms proportional to k 2 and |fc| 1+2a are 
still given by the expressions above (B.29). This yields again as above (B.29) 

b a 



F a (k) = 1 - "fk 2 + c a k 1+2a + 0(k 5 ) . 



(B.30) 



Appendix C. Explicit expression of F a (Y) for a = 1 

In this appendix, we give an explicit expression of Fi(Y) for a — 1. The starting point of our analysis 
is the expression ( [45] ) given in the text: 

' 2 (CI) 



Fi(Y) = 
1 



7T 1 + 4Y 2 

2(1 -8Y 2 ) 4 
7T V(l + 4F 2 )(l + 16y 2 ) + (l + 16y 2 )t 
For z a complex number, the following elementary relations are useful : 



Re 



(1 - AiY)i arctan ((1 + 4iY)~i 



arctan z 



1 

2~i 



(log(l + w)-log(l-«)) , 



log (x + iy) — log (\/x 2 + y 2 ) + 2i arctan [ J 

\x+^x 2 + y 2 



On the other hand has 
1 



VI + 4iF 



a + ib , a 



(i + ler 2 ) 1 ^ 
-i 



cos, - 



Defining A and /x as 

A = arctan 



(l + ier 2 ) 1 ^ 8111 v 2 

9 = arctan (4Y) . 
+ b 2 - 1 + v /a 2 + (l-&) 2 v /a 2 + (l + 6) 2 



2a 



/' 



- log 
4 B 



'{l-b) 2 +a 
{l + b) 2 +a 2 
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(C.2) 
(C.3) 

(C.4) 



(C.5) 
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in terms of a, b denned above (C.4), one obtains finally (after straightforward algebra) 

1 2 
7T f + AY 2 

2(1 -8Y 2 ) 4 / (W 

I A ens I — 



Fi(Y) = 



(C.6) 



1 



7T V(l + 4F 2 )(1 + 16F 2 ) (1 + 16F 2 ) 3 / 4 



/ism 
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